AD/A-004  226 


ACOUSTIC  FLUCTUATIONS  DUE  TO  SHALLOW 
WATER  THERMAL  MICROSTRUCTURE 

James  Burgess  Hagen 

Naval  Postgraduate  School 
Monterey,  California 

September  1974 


r 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  or  THIS  RAPE  (Whmi  Pit  I nlotod) 

REPORT  DOCUMENTATION  PAGE 


portnumI 


i  AftE  '  READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 
!.  00 VT  ACCESSION  NO.  A  RECIPIENT'!  CATALOO  NUNOCN 


«.  TITLE  (m  l  IvHItlo) 

Acoustic  Fluctuations  Due  to  Shallow 
Water  Thermal  Microstructure 


R.  TYNE  ON  REPORT  A  RERIOO  COVERED 

Master's  Thesis; 
September  1974 _ 

A  FERFORMINOORQ.  REPORT  NUMEER 


authors; 


James  Burgess  Hagen 


JMTRACT  OR  ORANT  NUMSERfe) 


».  PERFORMING  ORGANIZATION  NAME  AND  AOOREIS 

Naval  Postgraduate  School 
Monterey,  California  939^0 


JONAH  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMRERS 


II.  CONTNOLUNQ  OFFICE  NAME  AND  ADORERS 

Naval  Postgraduate  School 
Monterey,  California  939*10 


II.  REPORT  DATE 

September  197*» 

IS.  NUHRER  OF  PAOES 

jjrtr  15 2 


14.  MONITORINO  AOENCY  NAME  A  ADDRESS?!/  dllhrml  horn  CmMIIA*  OMImJ  IS.  SECURITY  CLASS.  (ol  Ml*  r*p*rtj 

Naval  Postgraduate  School  Unclassified 

Monterey,  California  9 39^0 

U*.  declassification/ OOWNORAOINA 
schedule 


IS.  OISTRIRUTION  STATEMENT  (ol  Ml*  R*p*rtJ 


Approved  for  public  release;  distribution  unlimited. 


IT.  OISTRIRUTION  STATEMENT  (ol  Iho  tkotraei  m  lond  In  Block  30.  II  dlllotonl  I 

Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 


US  Department  of  Commerce 
Springfield,  VA.  221SI 


I  IS.  supplementary  notes 


rngnmim  id  am 


IS.  KEY  WOROS  (Continue  on  teeefie  tide  II  nocotooiy  ltd  Idonilty  ky  Mock  numkoe) 

Acoustic  Fluctuations 

Thermal  Microstructure 

Shallow  Water  Thermal  Microstructure 


I  to.  A9STRACT  (Cfiilmf  m  §I4§  If  n«c«M«r  m4  Utility  ky  kl9  clr  numbw) 


Simultaneous  measurements  of  waves,  orthogonal  water 
particle  velocities,  temperature,  conductivity  and  acoustic 
phase  and  amplitude  measurements  at  20  and  *10  kHz  over 
horizontal  and  vertical,  2m  paths  were  made  at  the  NUC 
oceanographic  tower,  San  Diego  in  April  197*1. 

Four  stochastic  models  for  acoustic  propagation  including 
geometric,  Born,  Rytov  and  Debye  approximation  models  were 


nn  ro"M 

UU  I  JAN  71 

(Page  1) 


EDITION  OF  I  NOV  SS  IS  ORSOLCTE 
S/N  0 102*014*  (60 1  I 


UNCLASSIFIED _ 

I  SECURITY  CLASSIFICATION  OF  THIS  PAOE  fWNpo  DM 


UNCLASSIFIED 


-^JT9fn^ IW"*^1^  »P ’1  C^irr? —»•«*.  wywrw'.f  ■ 

i 


ftCUWITV  CL»IWFIC*TIOM  or  THU  POOtnriww  D.<« 


(20.  ABSTRACT  continued) 


tested.  The  Krosil'nikov  -  Obhukov  formulation  of  the  Born 
approximation  yielded  reasonable  results.  Dlffulty  In 
establishing  scales  obviated  food  results  from  the  Rytov  and 
Debye  models. 

Taylor  thermal  microscales  from  spatial  correlation 
analyses  determined  vertically  across  the  thermocllne  varied 
from  6.8  to  178  cm.  The  measured  variance  of  the  Index  of 
refraction  ranged  from  4.5X10”11  to  2.7xl0"3. 

Temporal  stochastic  analyses  showed  horizontal  acoustic 
path  losses  were  ten  times  vertical  losses  because  of 
diffraction  effects.  The  primary  mechanisms  for  acoustic 
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ABSTRACT 


Simultaneous  measurements  of  waves,  orthogonal  water 
particle  velocities,  temperature,  conductivity  and  acoustic 
phase  and  amplitude  measurements  at  20  and  40  kHz  over 
horizontal  and  vertical,  '2m  paths  were  made  at  the  NUC 
oceanographic  tower,  San  Diego  in  April  1974. 

Four  stochastic  models  for  acoustic  propagation  including 
geometric,  Bom,  Rytov  and  Debye  approximation  models  were 
tested.  The  Krosil^ikov  -  Obhukov  formulation  of  the  Bom 
approximation  yielded  reasonable  results.  Difficulty  in 
establishing  scales  obviated  good  results  from  the  Rytov 
and  Debye  models. 

Taylor  thermal  microscales  from  spatial  correlation 
analyses  determined  vertically  across  the  thermocline  varied 
from  6.8  to  1?8  cm.  The  measured  variance  of  the  index  of 
refraction  ranged  from  4.5x10  to  2.7x10  . 

Temporal  stochastic  analyses  showed  horizontal  acoustic 
path  losses  were  ten  times  vertical  losses  because  of 
diffraction  effects.  The  primary  mechanisms  for  acoustic 
variability  were  internal  and  surface  waves  pumping  the 
thermocline  up  and  down. 
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I.  INTRODUCTION 


A.  HISTORY 

A  continuing  program  of  In  altu  ocean  ml c restructure 
research  has  been  conducted  at  the  United  States  Naval 
Postgraduate  School  (NPS)  since  1971.  Under  this  program, 
a  series  of  three  annual  experiments  dedicated  to  the  study 
of  the  complex  temporal  and  spatial  relationships  between 
multiple  parameters  measured  In  shallow  water  have  been 
conducted  by  the  Department  of  Oceanography.  Each  of  these 
experiments  has  been  restricted  In  scope  to  measurements 
Inside  a  9  m^  volume  of  water  and  to  oceanographic  phenomena 
which  occur  In  the  frequency  range  0.001  to  2.5  Hz. 

B.  OBJECTIVES 

In  this  thesis.  Initial  analysis  of  certain  parameters 
measured  In  the  latest  of  these  experiments  Is  reported. 
Primary  emphasis  Is  given  to  stochastic  temperature  relation¬ 
ships  affecting  high  frequency  (20  kHz)  CW  sound  transmission 
over  a  short  path.  The  particular  objectives  of  this 
experiment  were  to  investigate  the  dependence  of  acoustic 
phase  and  amplitude  fluctuations  on  the  magnitude  and 
characteristics  of  the  ocean  thermal  structure  and  the 
Intensity  of  the  related  turbulence. 


II.  THEORY 


The  description  of  oceanic  turbulence  by  stochastic 
techniques  has  followed  developments  In  the  atmosphere. 

One  of  the  principal  difficulties  In  application  of  the 
theoretical  methods  provided  from  atmosj|Bric  research  has 
been  determination  of  the  range  of  applicability  of  certain 
turbulence  models  In  terms  of  ocean  parameters.  Another 
has  been  the  application  of  these  models  to  sound  propagation 
In  a  manner  similar  to  the  relationships  developed  for 
electromagnetic  wave  propagation  in  the  atmosphere  as 
summarized  by  Strohbehn  [1968].  Historically,  attempts  tc 
resolve  these  difficulties  have  been  interconnected  since 
their  primary  relevance  is  to  acoustic  prediction  which 
motivated  the  research. 

A.  STOCHASTIC  MODELS  FOR  THE  MEDIUM 

* 

Three  expected  value  techniques  which  have  been  applied 
to  turbulence  calculations  in  the  ocean  are  1)  the  spatial 
correlation  (autocorrelation)  function,  2)  the  structure 
function  and  3)  a  "piecewise  correlation"  function. 
Application  of  each  technique  requires  a  statistically 
homogeneous  medium  but  the  order  above  indicates  their 
relative  applicability  to  increasingly  less  stringent 
stationarity  requirements.  The  function  form  for  each  is 
listed  below. 
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The  autocovariance  function,  Rx(p): 

R x(p)  -  E[T(x)T(x+p)J.  .  (1) 

The  structure  function,  Dx(p) 

Dx(p)  »  E{[T(x)  -  Ttx+p)]2)  .  (2) 

The  piecewise  correlation  function,  R^(p) 

„r,  E{[T(x)  -  TCx+p' ) ]2} 

Ri(p)  E[1 - i  yo) 

D.,(p') 

" E[1  •  s"0'oy]  •  (3) 

A 

E  indicates  expected  value  and  p  is  the  spatial  lag. 

For  completely  stationary  processes,  the  structure  function 
and  the  "piecewise  correlation"  function  reduce  to 

Dx(p)  "  2ERX(0)  ~  Rx(p)]  »  (4) 

and  R£(p)  =  E[1  -  Dx(p) */2Rx(0) ]  -  Rx(p» )/Rx(0) . 

(5) 

The  stationarity  difficulties  often  encountered  in  the 

use  of  R  (p)  are  partly  overcome  by  the  use  of  D  (p)  because 
x  * 

the  large  scale  effects  in  the  ocean  are  differentially 
filtered  out  [Black  1965].  The  piecewise  correlation  function 
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can  be  applied  to  measurements  which  have  only  "piecewise 
atationarity."  This  makes  it  the  most  powerful  tool  of  the 
three  techniques.  "Piecewise  stationarity"  is  defined  as 
atationarity  (in  the  ordinary  sense)  over  a  sufficient 
record  length  for  which  the  structure  function  and  variance 
can  be  calculated.  This  method  was  proposed  by  Whitmarsh, 
Skudrzyk  and  Urick  [1957]. 

Various  analytical  equations  have  been  used  to  approxi¬ 
mate  the  form  of  the  auto  covariance  function  used  to  describe 
the  environment.  Llebermann  [1951]  used  a  simple  exponential 
form  in  his  early  experiments.  Others  have  used  a  Gaussian 
form  with  some  success  [Stone  and  Mintzer  1953].  Use  of  the 
exponential  form  yields  the  best  fit  to  experimental  data 
above  certain  spatial  lags,  but  it  underpredicts  the  number 
of  "small"  turbulent  blobs.  This  undesirable  effect  is 
eliminated  by  the  use  of  the  Gaussian  form. 

Kolmogorov  [19*11]  predicted  a  2/3  power  law  dependence 
for  the  structure  function  for  turbulent  flow  in  which  there 
was  a  decoupling  of  the  energy  producing  and  energy  diss..pa~ 
tion  regions.  Batchelor  [1953]  transformed  this  result  and 
demonstrated  that  in  wave  number  space  energy  spectra  have 
a  wave  number  to  the  -5/3  power  law  dependence.  This  was 
first  verified  convincingly  by  Grant  et  al.  [1962]  for  ocean 
measurements  and  Pond  et  al.  [1963]  for  atmospheric  measure¬ 
ments.  Ocean  measurements  do  not  always  show  the  dependence 
since  the  ocean  tends  to  be  stratified  much  more  intricately 
structured  than  the  atmosphere  which  makes  the  recording  of 
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continuous  stationary  data  unlikely.  However  the  -5/3 
power  law  depends  not  so  much  upon  isotropy  as  upon  the 
clear  separation  of  scales  of  energy  production  and  energy 
dissipation  In  which  a  cascade  of  energy  from  larger  to 
smaller  scales  occurs.  Whether  or  not  such  a  form  Is 
ubiquitous  In  the  ocean  is  not  as  significant  as  the  presence 
of  horizontal  sheets  which  could  cause  multiple  reflection 
layers.  Very  sharp,  multiple  thermal  layers  have  been 
reported  in  the  Arctic  [Denner  1971],  for  example. 

Recently  it  has  been  shown  that  surface  wave  effects  can 
be  removed  from  ocean  turbulence  measurements  with  the  result 
that  purely  turbulent  flow  remains  and  a  -5/3  power  relation¬ 
ship  may  be  observed  [Thornton  and  Boston  107*0.  These 
findings  were  based  upon  measurements  both  in  and  around  the 
thermocline  and  showed  that  anisotropy  was  not  a  critical 
difficulty,  provided  additional  energy  sinks  or  sources  did 
not  result  (e.g.  in  the  thermocline). 

Unfortunately,  other  energy  contributions  to  the 
turbulence  spectra  do  exist.  Microstructure  related 
turbulence  from  salt-fingering  is  one  form  [Gregg  1973, 

Denner  1971].  Billow  turbulence  resulting  from  shear 
instabilities  related  to  internal  waves  has  also  been 
reported  by  a  number  of  investigators  [Woods  and  Wiley  1972, 
Gregg  1973,  Powell  197*0*  Further,  the  existence  of  the 
internal  wave  itself  could  be  expected  to  add  energy  at  low 
frequencies  as  could  surface  wave  interaction  with  the 


thermocline. 


B.  WAVE  PROPAGATION  MODELS 

Pour  methods  derived  from  the  Helmholtz  equation  have 
been  used  to  describe  wave  propagation  through  the  turbulent 
ocean.  These  are:  1)  geometrical  optics,  2)  the  Bom 
approximation  [Bom  1933] »  3)  Rytov’s  method  [Rytov  19373, 
and  Jl)  the  Bebye  approximation  [Frisch  1968],  While  it  is 
not  a  necessary  assumption,  most  applications  of  these 
methods  presume  that  the  ocean  is  homogeneous  so  that  the 
scalar  wave  equation  can  be  used.  Such  an  assumption  places 
a  rather  severe  restriction  on  the  acoustic  wavelength,  X, 
such  that 


One  of  the  first  experiments  which  used  a  stochastic 
technique  to  examine  the  acoustical  effects  of  small-scale 
temperature  lnhomogeneltles  In  the  ocean  was  conducted  by 
Llebermann  [1951].  Through  the  use  of  the  classical  optical 
model,  Llebermann  was  able  to  predict  certain  scintillations 
In  sound  amplitude  observed  at  sea  by  the  equation: 


Vj^  1  (L/a)^ 


Vj2  ■  E[(^)2] 


E[<^)2] 

co 


where  Vj  ■  variance  of  the  sound  intensity,  IQ  ■  mean  sound 
intensity,  AI  ■  the  instantaneous  deviation  of  sound  Intensity 

p 

from  the  mean,  o  ■  variance  of  the  index  of  refraction, 

L  *  acoustic  path  length  and  a  ■  Taylor  temperature  micro¬ 
scale  [Lumley  and  Panofsky  196*1  ].  Llebermann  also  showed 
that  sound  reverberation  encountered  in  the  ocean  was  larger 
than  that  predicted  due  to  thermal  inhomogeneities  alone. 

Thid  was  substantiated  by  later  developments  which  have 
shown  that  biological  scattering,  chemical  relaxation  and 
bubbles  in  the  near-surface  region  are  important  contributors 
to  sound  attenuation. 


Sagar  [1955,1957,1959]  was  able  to  show  that  sonar 
amplitude  fluctuations  resulting  from  temperature  measure¬ 
ments  could  be  correlated  to  sea  state.  Additionally,  he 


showed  that  the  linear  dependence  of  the  variance  of  sound 


pressure  level  on  range  [Sheehy  1950]  was  apparent  only  when 
thermal  microstructure  was  present. 

A  modified  geometrical  optics  technique  has  recently 
been  proposed  by  Tatarskl  [1967].  The  method  Is  similar  to 
the  classical  one  except  the  path  length  Is  described  as 
the  elkonal  or  total  path  length,  8.  However,  one  would 
assume  that  this  would  modify  the  restriction  on  the  results 
to 


a  «  *o2/e  (ii) 

No  application  of  this  method  to  the  ocean  has  yet  been  made. 
2.  Born  Approximation 

The  Born  approximation  is  a  perturbation  technique 
applied  to  the  wave  equation  which  is  valid  only  as  long 
as  the  acoustic  amplitude  fluctuations  remain  small.  This 
is  not  a  serious  restriction  in  the  ocean,  even  under  the 
most  severe  scattering  conditions,  but  the  model  only  accounts 
for  a  single  scattering  from  a  turbulent  blob  over  the 
acoustic  path  length,  L.  This  requires  that 

A  <<  o  /aL  (12) 

The  geometrical  optics  results  are  a  special,  high  frequency 
case  included  in  this  model. 


Stone  and  Mlntzer  [1962,  1965]  conducted  a  series 
of  laboratory  experiments  from  which  they  demonstrated 
solutions  to  the  Born  Model  which  predict  the  variance  of 
normalized  sound  pressure  level,  V2,  in  a  iow  frequency 
region  (wave  region)  and  in  the  ray  theory  region.  For 
the  wave  region  they  showed: 

V2  -  1  kQ2  o2  a  L  (Gaussian  Rx(p)  (13) 

V2  -  a2  kQ2  a  L  (Exponential  Rx(p)  (1*0 


V2  -  E[(^)2]  (15) 

o 

X  «  2ira  (16) 

X  «  2irL  (17) 

kQ  -  2ir/X  (18) 


Where  PQ  ■  mean  sound  pressure  level  and  AP  ■  pressure  level 
difference  from  the  mean.  For  the  ray  region,  no  solution 
is  possible  if  the  exponential  form  of  the  spatial  correlation 
function  is  used  since  the  solution  is  divergent.  But  they 
showed  that  the  f<  ilowing  equation  results  from  use  of  the 
Gaussian  correlation  function: 
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A  a^(Va)^ 


(19) 


Neubert  and  Lumley  (1970]  used  data  from  the  above 
experiment,  to  eomp.re  with  atmospheric  data  baaed  upon  the 
..w  model.  Their  result,  showed  that  the  turbulent  scale, 
for  air  and  water  are  of  the  same  order  of  magnitude  for 
"similar"  turbulent  flow..  They  were  also  able  to  better 
•  define  the  limit  on  wave  length  specified  in  relation  W, 

such  that 

X<a/5[ok(u/c0)]? 
a2  -  Xg2/ok 

Rt  * 


(21) 

(22) 


where : 

X  *  Taylor  turbulence  microscale, 
g 

o  ■  Prandtl  number  for  the  medium, 
k 

u  ■  local  turbulent  velocity, 
c  -  mean  speed  of  sound  in  the  medium, 


Richardson  number, 


and 


n  ■  Kolmogorov  turbulence  microscale. 

8 

In  addition,  they  found  a  restriction  on  acoustic  frequency 
was  necessary,  though  it  was  usually  less  restrictive  than 
(20).  This  restriction  is 

1  1 

f  >  750[ok(Rt)7]?/Xg  (sec"1)  (23) 

It  follows  from  (6)  and  (20)  that 

1 

lQ  z  0{a/50[ok(u/co)]?}  (cm)  (24) 

which  lends  additional  meaning  to  relations  (6),  (7)  and 

(11). 

Whitmarsh,  Skudrzyk  and  Urick  [1957]  conducted 
experiments  in  the  ocean  which  largely  agree  with  Stone’s 
and  Mintzer’s  laboratory  results.  However,  they  noted  the 
variance  in  sound  pressure  level  due  to  thermal  scattering 
varied  directly  as  the  acoustic  frequency  and  as  the  range 
squared  in  the  Kolmogorov  turbulence  range.  This  result  is 
obviously  in  the  intermediate  range  between  the  ray  and 
wave  solution  regions  suggested  by  Stone  and  Mintzer. 

They  also  observed  that  the  turbulent  patches  were 

s 

/ 

not  spherical,  but  were  spheroidal  in  shape.  In  this  case, 
a  simple  transform  to  cylindrical  coordinates  shows  that  the 

\ 
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Ttylor  microscale,  a,  must  be  replaced  such  that 

a2  -  (b2/3)(2  +  c2/b2)  (25) 

where  b  ■  narrower  dimension  and  c  -  longer  dimension  of 
the  turbulent  patch.  This  has  the  effect  of  reducing  the 
focusing  length  of  an  individual  turbulent  "lens"  which,  ‘ 
in  turn,  increases  the  effective  ray  solution  region. 

The  following  expression  for  the  variance  of  the 
sound  pressure  level,  V2,  has  been  presented  as  valid  over 
the  entire  wave  and  ray  regions  [Krasil'nikov  and  Obhukov, 

1956[: 

V2  -  kQ2  o2  a  L[1  -  l/Dtan"1D]  (26) 

D  ■  irLAQa.  ^2^ 


For  this  result,  the  wave  region  is  defined  when 


X  »  2a2/L 


(28) 
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3.  Rytov^  Method 

Rytov’s  method  (method’ of  smooth  perturbations)  is 
the  method  used  by  Tatarskl  [1961]  to  develop  a  stochastic 
wave  model  for  the  atmosphere.  It  includes  the  geometrical 
optics  and  Bom  approximation  methods  as  special  cases  and 
includes  multiple  scattering  effects  as  well.  Unfortunately, 
the  scattering  effects  are  Included  such  that  they  cannot 
be  parameterized  directly  in  terms  of  the  turbulence. 

Either  the  spatial  correlation  or  structure  function 
may  be  used  with  this  method  as  dictated  by  considerations 
regarding  the  statlonarlty  of  measurements  to  be  analyzed. 

The  restriction 

X  «  YLo/L  (30) 

is  necessary  for  reasonable  solutions  of  the  resultant 
equations.  It  Is  also  necessary  to  compute  the  range  of 
validity  of  the  perturbation  expansion  when  this  method  Is 
used.  This  calculation  will  not  be  discussed,  but  Strohbehn 
[1968]  suggested  that  the  method  Is  valid  as  long  as  the 
variance  In  the  normalized  log-amplitude  of  the  sound  pressure 
level  Is  less  than  0.8.  He  also  argued  that  the  results 
from  sound  phase  calculations  should  be  less  restricted  than 
those  for  amplitude. 

Black  [1965]  applied  this  model  to  certain  data  from 
the  CHAIN  experiment.  One  inportant  conclusion  he  reached 
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was  that  ths  standard  deviation  of  the  fluctuation  in  sound 


Intensity  through  a  turbulent  medium,  aa,  and  the  principal 
period,  Tf,  of  its  fluctuation  can  be  predicted  by  the  equations: 


17  -  A 

0#  -  12.7  Cn  L12  X  12 


Tb  •  ir/n:/2vn  (32) 

C2n  -  Ct2/HT02  +  Cs2/1000S2  +  Cy2/c02  a  Ct2/4T0  (33) 

2 

CT2  -  hj(L)/L^  (33a) 

2 

CB2  -  D8(L)/^  (33b) 

2 

Cv2  -  DV(L)/L*  (33c) 


where : 

v  *  absolute  value  of  the  medium’s  mean  velocity 
across  the  acoustic  path, 

TQ  *  mean  absolute  temperature  (°C), 

% 

S  ■  mean  salinity  in  parts  per  thousand, 

CQ  ■  mean  speed  sound  (m/sec), 

Dl(L)  ■  value  of  the  temperature  structure  function 
at  range  L, 

D  (L)  ■  value  of  the  salinity  structure  function 
at  range  L, 

D  (L)  ■  value  of  the  turbulence  structure  function 
v  at  range  L. 


The  values  for  C„  which  Black  observed  were: 
n 

_ij  _i 

In  thermocline:  5.2  x  10  m 

-4  -i 

Hear  surface:  1.0  x  10  m 

_C  _1 

Below  thermocline:  2.1  x  10  'm  . 
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4*  Debye  Approximation 

Tht  Debye  approximation  to  the  wave  equation  la  the 
■oat  general  model  currently  available  for  the  deaorlptlon 
of  the  turbulent  ocean  effecta  on  eound  tranamiaalon  [Neubert 
1970].  It  contalna  Rytov'a  method  aa  part  of  it a  aolutlon 
and  It  reducea  to  the  Bom  model  when  the  alngle  acatterlng 
condition 

X  »  2  1  a/in  (34) 

la  applied.  The  reaultant  equation  from  the  model, 

V2  «  2  a2  kQ2  Le  L  (35) 

where  L_  11  the  Eulerlan  Integral  scale  (e.g.  the  largest 
turbulent  scaled  has  no  special  long  range  nor  high  frequency 
restrictions.  In  other  words,  it  is  good  when  the  cumulative 
effects  of  small,  continuous  fluctuations  of  the  local 
turbulent  velocity  strongly  affect  phase  at  longer  ranges. 

No  use  has  been  made  of  this  model  in  the  ocean  to  date. 

C.  MODEL  APPLICATION 

Stewart  and  Grant  [1962]  conducted  a  turbulence  experiment 
near  the  surface  in  the  presence  of  large  waves.  They  showed 
the  Taylor  microscale,  a,  was  on  the  order  of  6  cm  which 
corresponds  to  the  results  (6.4  cm)  Neubert  and  Lumley  [1970] 
calculated  from  Sheehy's  [1950]  data.  Since  these  measure¬ 
ments  correspond  well  and  fit  the  conditions  of  the  current 
NPS  experiment,  the  Neubert  and  Lumley  values  were  used  to 
provide  a  rough  basis  for  further  analysis.  These  values, 
coupled  with  experimental  parameters,  are  shown  in  Table  1. 


The  limits  of  applicability  of  each  of  the  wave  models  and 
the  associated  value  of  V2  calculated  on  the  basis  of  this 
data  Is  shown  In  Table  2.  It  would  appear  from  the  latter 
table  that  the  Debye  approximation  and  perhaps  the  Bom 
model  for  the  wave  region  (If  the  limit  Imposed  by  relation 
(6)  Is  not  too  strong)  would  be  most  applicable  In  this 
ease.  The  large  difference  In  V2  resultant  from  the  Bom 
model  compared  to  that  from  the  Debye  model  could  be  the 

t 

result  of  poor  definition  of  the  micro  and  Integral  scales. 

One  should  also  note  that  the  acoustic  near-field 
criteria  for  a  piston  source  Imposes  an  additional  requirement 
on  the  wavelength: 

X  >  r2/L  (36) 

where  r  *  radius  of  the  transmitter  active  disc.  Since 
relation  (36)  requires  X  >  0.6  cm,  geometrical  optics  and 
Rytov's  method  can  be  immediately  discarded. 


TABLE  1 

PARAMETER  VALUES  FOR  MODEL  TESTING 


Parameter 

Assumed  Value 

Source 

c 

1505  m/sec 

Experiment 

0 

0 

0 

7  x  10”5 

NL 

L 

40m 

NL 

e 

A 

6.4  cm 

NL 

g 

6300 

NL 

°k 

7 

NL 

«2a 

5.0  x  10”7 

S 

L 

200  cm 

Experiment 

f 

20  khz 

Experiment 

A 

7.5  cm 

Experiment 

c 

-4  -1 

1.0  x  10  m  A 

B 

n 

u 

10.5  cm/sec 

C 

nE 

258  cm 

Equation  22 

a 

2.4  cm 

Equation  21 

*0 

2.2  cm 

Equation  24 

a 

4.6  x  10"4 

C 

It 

0.84  cm”'1' 

Equation  18 

0 

D 

130 

Equation  27 

r 

10.88  cm 

Experiment 

NL  *  Neubert  and  Lumley  [1970] 

S  ■  Sheehy  [1950] 

B  ■  Black  [1965] 

C  ■  Calculated  from  above  data. 
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TABLE  2 


RESULTS  FROM  MODEL  TESTING 


Model 

ifeat  Restrictive  Limit 

Geometric 

5.8  x  10"2  (19) 

X  «  1.1  x  10"2cm  (7) 

Bom 

(Gaussian  fora) 

6.3  x  10"5  (13) 

1.0  x  10~2 «X«  15cm(12,l6) 

Bom 

(Exponential  form) 

7.1  x  10-5  (14) 

1.0  x  10"2  «  X  «  15cm(12,l6) 

Bom 

(Neubert  Limits) 

X  <  20.9cm  (18) 

Bom  (K— 0) 

7.1  x  10"5  (26) 

X  »  2.9  x  10“2cm  (28) 

Rytov 

1.1  x  10“2  (31) 

X  «  0.5  cm  (30) 

Debye 

(L  ■  40m) 

© 

2.4  x  10"1  (35) 

None 

Debye 
(L  -  3  m) 

1.8  x  10'2  (35) 

None 

HOMOGENEOUS  ASSUMPTION 
X  «  2.2cm  (6) 

AflOIISTIC  FAR-FIELD  CRITERION 
X  >  0.6  cm  (36) 

Note:  Numbers  in  parenthesis  refer  to  applicable  equation. 
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III.  EXPERIMENT  AND  DATA  ANALYSIS 


A.  EXPERIMENT 

Measurements  were  taken  from  May  24  to  26,  197^  from 
the  Naval  Undersea  Research  and  Development  Center’s  (NUC) 
oceanographic  research  tower  located  1.6  km  off  Mission 
Beach  San  Diego,  California.  The  tower  is  a  stable,  three¬ 
tiered  concrete  and  steel  structure  built  into  sandy  bottom 
in  water  of  18.5  m  mean  depth.  It  provides  enclosed  spaces 
for  electronic  equipment  and  open-deck  working  areas  for 
handling  instruments.  Winch  controlled  trolleys  are  avail¬ 
able  on  three  sides  of  the  structure  on  which  devices  may 
be  lowered  to  within  1  m  of  the  bottom.  Instruments  for 
this  experiment  were  installed  on  an  aluminum  and  steel 
frame  attached  to  a  trolley  on  the  west  side  of  the  tower 
which  faced  the  predominant  swell.  A  schematic  view  of 
the  experiment  layout  is  shown  in  Figure  1. 

Basic  meteorological  and  oceanographic  data  were  avail¬ 
able  from  equipment  permanently  installed  on  the  tower.  A 
bathythermograph  trace  obtainable  from  the  north  side  of 
the  structure  facilitated  placement  of  the  instrument  frame 
in  the  desired  vertical  position  with  respect  to  the 
thermocline. 

The  general  environmental  conditions  encountered  were 
moderate  winds,  sea  and  swell  from  the  West.  Light  winds 
to  5  kts,  overcast  skies  and  swell  to  1  m  were  usual  in  the 


■oming  (set  Table  3)*  In  the  afternoon,  winds  increased 
to  15  kts  and  seas  of  0.75  ra  with  accompanying  swell  of  1  m 
was  common.  The  diurnal  variation  In  thermocllne  depth  was 
approximately  12  m. 

1.  Instruments 

Temperature  measurements  were  made  with  Fenwell 
Electronics  Inc.  type  k496,  double-bead  thermistors  which 
exhibit  low  self-heating  and  well-matched,  high-resolutlon 
characteristics.  Calibration  of  the  thermistors,  together 
with  associated  circuitry,  was  conducted  In  the  laboratory 
under  closely  controlled,  stable  conditions.  A  Hewlett- 
Packard  HP-2801A  quartz  thermometer  with  an  accuracy  of 
±  ,006°C  was  used  as  a  temperature  standard  and  a  decade 
box  was  used  as  a  resistance  standard.  The  same  decade  box 
Was  used  to  calibrate  the  thermistor  array  in  the  field  to 
eliminate  resolution  difficulties  resulting  from  the  poor 
overall  accuracy  of  the  decade  box. 

Water  particle  velocities  were  measured  using  Marsh- 
McBirney  electromagnetic  flowmeters  (models  721  and  722). 
Since  these  flowmeters  were  calibrated  by  the  manufacturer 
only  for  steady  flow  conditions,  each  probe  was  recalibrated 
for  unsteady  flow  by  oscillating  it  in  a  water  tank  at  rates 
corresponding  to  expected  experiment  velocities. 

A  Baylor  Co.  wave  staff  system,  model  13528R,  was 
used  to  measure  sea  surface  elevation.  Calibration  was 
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accomplished  at  the  NUC  tower  by  successively  shorting  out 
sections  of  the  sensor  active  element  while  measuring  output 
voltage . 

A  Naval  Research  Laboratory,  Underwater  Sound 
Reference  Division  (USRD)  type  P-33  general  purpose  uni¬ 
directional  transducer  was  used  as  a  CW  sound  transmitter. 

It  features  two  piezoelectric  ceramic  arrays  mounted  coaxially 
which,  when  excited  in  parallel,  exhibit  exceptional  direc¬ 
tivity.  This  was  desired  to  reduce  interference  from  surface, 
bottom  and  equipment  frame  reflections.  The  beam  widths 
(— 3db  down)  for  20kHz  and  40kHz  were  20°  and  10°  respectively. 
USRD  type  P-50  lead  zirconate-titanate  omnidirectional 
transducers  were  used  as  receivers. 

The  sonar  array  was  calibrated  in  an  anechoic  test 
tank  filled  with  fresh  water  and  the  hydrophone  calibrations 
compared  satisfactorily  with  calibrations  previously  per¬ 
formed  at  USRD.  .  The  frequency  standard  used  was  a  Hewlett- 
Packard  HP-2801A  quartz  thermometer  configured  as  a  frequency 
meter  and  the  voltage  standards  were  an  RMS  calibrated  Fluke 
digital  voltmeter  and  a  Hewlett-Packard  model  HP  3400A  RMS 
voltmeter.  Phase  measurements  were  conducted  using  a  Dranetz 
model  305  phase  meter.  Calibration  and  accuracy  measurements 
for  the  peripheral  equipment  associated  with  the  sound  arrays 
were  taken  concurrently  with  transducer  calibration. 

The  instrumentation  accuracy,  response  and  limitations 
are  summarized  in  Table  4. 
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2.  Instrument  Arrangement 

Two  arrays  of  thermistors  were  mounted  on  one  side 
of  the  equipment  frame  In  a  common  vertical  plane.  Seven 
thermistors  were  mounted  along  a  horizontal  member  and  7 
were  mounted  along  a  vertical  member  as  shown  In  Figure  2. 
The  Instrumentation  and  abbreviated  notation  used  In  the 
text  Is  given  In  Table  5.  One  thermistor  In  the  vertical 
array  was  common  to  the  horizontal  array.  The  thermistor 
probes  were  spaced  In  an  approximate  geometric  progression 
so  that  a  maximum  number  of  spatial  combinations  could  be 
achieved  for  adequate  definition  of  the  spatial  correlation 
and  structure  functions  with  a  limited  number  of  probes. 

The  minimum  spacing  (5.0  cm)  was  determined  by  the  physical 
size  of  the  thermistor  mounts.  This  was  compatible  with 
the  smallest  wavelengths  which  could  be  measured  within  the 
time  constant  of  the  thermistor.  Maximum  thermistor  spacing 
was  determined  by  physical  constraints  on  the  length  of  the 
support  pipe  to  prevent  strumming. 

Three  velocimeters  were  placed  so  that  the  horizontal 
and  vertical  particle  velocities  could  be  measured  in  the 
same  plane  as  the  thermistor  arrays.  Additionally,  they 
were  placed  in  a  vertical  line  directly  under  the  Baylor 
wave  gage  with  the  horizontal  component  measured  from  the 
west.  The  temperature  and  particle  velocity  devices  were 
offset  29  cm  from  the  support  frame  by  mounting  brackets  to 
minimize  false  turbulence  effects  generated  by  the  frame. 
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FIGURE  2.  RELATIVE  POSITIONS  OF  INSTRUMENTS  ON  FRAME  (SDE  VIEW) 
INSERTS  SHOW  SONAR  ARRAY  CENTERLINE  POSITION  RELATP*’  TO  FRAME  BAVUOR 
WAVE  GAGE  WAS  LOCATED  DIRECTLY  ABOVE  VERTICAL  ARRAY  DEVICES. 


TABLE  5 

CHANNEL  ABBREVIATIONS 


n«t-a  Channel 

Abbreviation 

Order  in  Pig.  2 

y&Vll  VlW***®  * 

Vertical  Temperatures 

TMV1-TMV7 

From  Bottom 

Horizontal  Temperatures 

TMH1-TMH7 

From  West 

Sound  Amplitude 

SAM1-SAM3 

‘From  East 

or 

Sound  Phase 

SPH1-SPH2 

From  Top 

Sea  Surface  Elevation 

WAVE 

- 

Horizontal  Plow  Velocity 

PVAY-FVCY 

From  Bottom 

Vertical  Flow  Velocity 

PVAX-PVCX 

From  Bottom 

Conductivity 

SAL1-SAL2 

From  West 

The  sonar  array  was  separately  mounted  on  a  tripod 


base  which  pivoted  vertically  around  the  main  equipment 
frame  but  which  was  held  rigid  by  clamps  during  recording 
periods.  This  arrangement  facilitated  horizontal  and 
vertical  sound  path  measurement  (Figure  2).  The  centerline 
of  the  transmitted  sound  beam  was  57  cm  away  from  the 
vertical  plane  of  the  thermistor  arrays  to  minimize  unwanted 
reflections  from  the  other  Instruments. 

Three  F-50  receivers  were  placed  along  the  center 
of  the  transmitted  beam  outside  the  near  field  at  2.0,  2.3 
and  2.9  ra  from  the  transmitter  face.  The  minimum  separation 
(30  cm)  was  selected  to  maximize  the  number  of  wavelengths 
through  which  the  lowest  frequency  (20kHz)  sound  would 
travel  between  receivers  and  to  remain  near  the  lower  spatial 
scales  of  the  measured  temperature.  The  maximum  sound  path 
length  of  2.9  m  was  selected  to  minimize  reflections  from 
the  frame  and  mountings. 

3.  Pre-recording  Signal  Processing 

The  signals  from  the  wave  height  and  particle  velocity 
devices  were  recorded  directly.  All  other  signals  were 
processed  through  a  combination  of  filters,  amplifiers  and 
dc  bias  removal  devices  to  maximize  signal-to-noise  (S/N) 
ratio  before  recording. 

a.  Temperature  Circuitry 

Each  of  the  horizontal  and  vertical  temperature 
arrays  was  connected  to  a  separate  temperature  array  drawer 


which  contained  Wheatstone  bridge  networks,  amplifiers  and 
calibration  facilities  for  seven  channels  (see  Figure  3)* 

The  Wheatstone  bridge  network  provided  the  capability  to 
match  the  response  of  each  channel.  Additionally,  balance 
adjustments  were  provided  to  adjust  the  operating  point  of 
the  amplifiers  to  the  mean  temperature  so  that  distortion 
during  recording  was  minimized, 
b.  Sonar  Circuitry 

As  shown  in  the  block  diagram  (Figure  4),  signals 
from  the  F-50  transducers  were  amplified  and  then  band-pass 
filtered  at  the  operating  frequency.  The  signal  was  then 
processed  through  a  second-detector-type  demodulator.  The 
dc  bias  was  removed  from  the  detected  signal  and  the  ampli¬ 
tude  modulation  (AM)  signal  was  then  amplified  and  low-pass 
filtered. 

Dranetz  type  305  phase  meters  were  connected 
between  the  band-pass  filter  outputs  from  receiver  hydrophones 
one  and  two  and,  similarly,  between  hydrophones  two  and  three. 
The  dc  bias  was  removed  from  the  output  of  each  phase  meter 
and  the  signal  was  then  amplified  and  low-pass  filtered. 

Before  each  run  the  dc  bias  setting  was  made  and 
the  rms  voltage  from  each  band-pass  filter  output  and  mean 
phase  shift  read  from  each  phase  meter  was  recorded.  This 
technique  allowed  much  higher  signal  amplification  to  be 
used  prior  to  recording  because  the  signal  amplifiers  were 
not  saturated  by  the  large  dc  segment  of  the  signal. 
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FIGURE  l.  BASIC  TEMPERATURE  MEASUREMENT 

CIRCUIT  BLOCK  DIAGRAM.  ONLY  I  OF  THE  7  CHANNELS  IS 
SHOWN  IN  THE  TEMPERATURE  ARRAY  DRAWER  FOR  CLARITY. 
TWO  IDENTICAL  ARRAY  DRAWERS  WERE  USED  DURING  THE 

EXPERIMENT. 
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FIGURE  4.  BASIC  SONAR  PHASE  AND  AMPLITUDE 
MEASUREMENT  BLOCK  DIAGRAM.  TWO  PHASE  CHANNELS  AND 
THREE  AMPLITUDE  CHANNELS  WERE  USED,  HOWEVER  ONLY 
ONE  OF  EACH  IS  SHOWN  FOR  CLARITY. 
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4.  Experimental  Procedure 

Twelve  data  recording  runs  of  29  to  73  minutes 
duration  were  taken  at  various  depths  above,  In  and  below 
the  thermocline.  A  total  of  28  channels  of  data  were 
recorded  for  each  run  Including  2  channels  dedicated  tc 
experimental  conductivity  devices  [Kane  197*0.  Surface 
elevation,  vertical  and  horizontal  particle  velocities, 
temperature,  sound  amplitude  modulation,  sound  phase  modu¬ 
lation  and  conductivity  were  measured  simultaneously. 

Run  depths  were  selected  on  the  basis  of  a  BT  trace 
taken  just  before  the  run,  but  the  final  run  depth  was 
chosen  to  optimize  the  degree  of  fluctuation  in  temperatures 
monitored  from  the  vertical  thermistor  array.  The  sonar 
array  was  held  in  a  fixed  position  during  each  recording 
run.  Runs  were  made  with  the  sonar  array  oriented  hori¬ 
zontally,  vertically  and  at  an  angle  with  respect  to  the 
support  structure.  This  allowed  data  to  be  obtained  at 
various  attitudes  through  and  along  the  thermocline.  A 
summary  of  all  twelve  runs  is  given  in  Table  6. 

A  Sabre  III  and  a  Sangamo  model  3500,  both  1*1  channel 
FM  tape  recorders,  were  used  to  record  the  analog  data.  A 
timing  pulse  was  recorded  on  one  channel  of  each  recorder 
at  the  beginning  of  each  run  to  provide  temporal  continuity 
between  recorders. 

The  analog  recordings  were  digitized  using  a  hybrid 
system  including  a  Scientific  Data  System  model  5000  analog 


computer  and  a  Xerox  Data  Systems  model  9300  digital  computer. 
Digital  data  preparation  and  analysis  was  performed  on  an 
IBM  model  360  computer.  The  digitizing  data  rate  used  was 
25Hz,  however  60Hz  electronic  noise  was  coupled  into  the 
signal,  so  low-pass  filters  were  used  to  remove  it.  Addi¬ 
tionally,  noise  spikes  which  existed  in  some  channels  were 
removed  by  a  digital  editing  routine.  Data  handling  capacity 
was  improved  by  reducing  the  number  of  samples  by  five  using 
an  inverse  transform,  numerical  low-pass  filter  [Davidson 
1970].  The  latter  reduced  the  analysis  Nyquist  frequency 
to  2.5Hz.  The  linear  trend  and  remaining  dc  bias  were  also 
removed  before  analysis. 

Spectral  analysis  was  performed  on  26.6-  minute 
segments  of  certain  data  sets.  This  choice  dictated  a 
minimum  detectable  frequency  of  0.00625Hz.  A  standard  lag 
time  of  5%  of  the  record  length  was  chosen  which  resulted  in 
1|0  degrees  of  freedom  for  each  computed  spectrum.  The 
resultant  80?  confidence  limits  of  the  chi-squared  distri¬ 
bution  were  between  0.73  and  1.30  of  the  measured  spectral 
estimates.  A  Parzen  window  was  applied  to  the  correlation 
function  during  the  spectra  calculations  to  minimize  the 
effect  of  finite  record  length  on  the  spectral  analysis. 


IV.  EXPERIMENT  RESULTS 


During  this  experiment,  the  primary  factors  influencing 
the  results  were  the  sea  state,  the  thermocline  gradient  and 
the  position  of  the  instruments  with  respect  to  the  thermo¬ 
cline.  The  sea  state  remained  nearly  constant  throughout 
the  recording  periods  although  the  peak  of  the  swell  showed 
a  gradual  shift  to  higher  frequencies.  The  thermocline 
depth  and  gradient  continuously  fluctuated  under  the  influence 
of  internal  waves  which  are  characteristic  of  the  waters  off 
San  Diego  during  this  time  of  year.  In  order  to  resolve 
some  of  the  effects  of  the  thermolcine  variability,  analyses 
in  and  above  the  layer  were  required.  Experiment  runs  one 
and  two  were  chosen  for  this  purpose.  Additionally,  run 
12  was  examined  to  check  the  results  since:  1)  it  included 
data  in  and  above  the  thermocline  in  a  single  record;  2)  the 
thermocline  gradient  was  intermediate  between  the  gradient 
encountered  in  runs  one  and  two;  3)  the  sonar  array  was 
angled  toward  the  vertical  rather  than  horizontal  during  the 
run  and  4)  the  run  was  conducted  two  days  after  runs  one 
and  two  under  slightly  different  sea  conditions. 

A.  ENVIRONMENT 

The  thermocline  migrated  diurnally  from  near  the  surface 
in  the  morning  to  about  14m  in  the  evening.  It  was  generally 
diffuse  near  the  surface  but  near  the  bottom  it  was  compressed 
so  that  a  very  sharp  thermal  gradient  (2°C/m)  was  created. 


Superimposed  upon  this  effect  was  an  apparent  four-hourly 
dilation  and  contraction  of  the  Isotherms. 

Figures  5-9  show  representative  bathythermograph  (BT) 
traces  taken  periodically  during  the  experiment.  The  traces 
are  plotted  so  that  the  12°C  mark  corresponds  to  the  time 
the  BT  was  taken.  Recording  run  envelopes  In  depth  and  time 
are  included  to  show  the  position  of  the  Instrument  array 
relative  to  the  thermocline.  The  11°,  12°  and  13°C  isotherms 
are  shown  to  indicate  the  relative  intensity  of  the  tempera¬ 
ture  gradient  during  the  runs.  Figure  9  is  an  expanded 
time  plot  for  BT’s  taken  during  run  12  which  shows  a  large 
internal  wave  having  a  36  min  period;  this  was  observed 
throughout  the  experiment. 

The  four-hourly  variation  in  the  thermocline  gradient 
probably  was  caused  by  internal  wave  interference  effects 
since  analysis  showed  that  these  variations  were  associated 
with  periods  of  high  temperature  variance  but  not  with  sur¬ 
face  waves  nor  heat  budget  considerations. 

Surface  slicks  presumably  associated  with  internal  waves 
were  seen  propagating  from  the  West  at  certain  intervals 
during  the  experiment.  A  typical  example  of  the  effect  of 
these  waves  is  shown  in  the  vertical  temperature  analog 
records  of  run  number  12  (Fig.  10).  The  large  36  min  wave 
peak  passing  the  vertical  temperature  array  is  shown  by  the 
dramatic  temperature  decrease  which  occurred  at  about  1800. 
Nine  and  1.5  min  period,  small-amplitude  internal  waves 
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FIGURE  6.  BATHYTHERMOGRAPH  TRACES  FOR  THE  MORNING  OF  29  APRL 


FIGURE  7  .  BATHYTHERMOGRAPH  TRACES  FOR  THE  AFTERNOON  OF  25  APRIL. 


FIGURE  9  .  SHORT  TIME  INTERVAL  BATHYTHERMOGRAPH  TRACES  FOR  THE 
PERIOD  OF  RECORDING  RUN  12 ,  26  APRIL. 


effects  are  seen  superimposed  upon  the  larger  internal  wave. 
The  "grass"  is  surface  swell  related.  Spectral  analysis 
showed  that  a  broad  band  of  Internal  wave  periods  from 
160  secs  (maximum  detectable  by  the  method)  to  as  low  as 
30  secs  were  superimposed  on  the  large  wave  as  well. 

Figure  11  shows  the  resultant  sound  effects  from  the 
internal  wave  temperature  perturbations  of  Fig.  10  during 
run  12.  The  sound  array  was  angled  25°  to  the  vertical  in 
run  12.  The  sound  amplitude  modulated  signals  SAMI  and  SAM2 
were  initially  above  the  layer  and  were  not  affected  by 
the  turbulent  thermocline  losses  as  severely  as  SAM3  (below 
the  layer).  As  the  layer  moved  up,  the  interference  and 
diffraction  effects  of  the  turbulent  blobs  [Whitmarsh  and 
others  1957]  in  the  vicinity  of  SAM3  were  diminished.  This 
is  shown  by  the  increase  in  sound  pressure  level  (about 
10-3  ub)  on  the  SAM3  analog  record. 

When  the  layer  moved  into  the  vicinity  of  SAMI  and  SAM2, 
the  sound  pressure  level  at  these  transducers  was  reduced 

„  _ji  _ii 

by  about  6  x  10  ub  and  3  x  10  ub,  respectively.  The 
thickness  of  the  "grass"  on  these  plots  is  an  indication  of 
the  wave-induced  amplitude  scintillation  effects  (~2  to  5 
x  10  ub)  caused  by  pumping  the  turbulent  blobs  back  and 
forth  across  the  sound  path.  (Similar  results  occurred  when 
the  sound  path  was  horizontal,  but  the  scintillations  were 
then  caused  by  waves  pumping  the  thermocline  up  and  down.) 
Additionally,  the  movement  of  the  layer  (AT  s  .3°C)  between 
SAMI  and  SAM2  caused  sound  phase  modulation,  SPH1,  to  phase 
shift  about  30°. 
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b.  tem:  oral  stochastic  analysis 


Tine  aeries  analyses  were  performed  on  26.67  minute 
segments  of  each  of  the  sound  channels,  the  wave  record, 
three  temperature  records  and  two  water  particle  velocity 
records.  Analyses  included  for  each  signal  the  calculation 
of  the  first  and  second  moments,  variance,  mean  period 
(calculated  from  total  zero  upcrossings) ,  autocorrelation 
function  and  autospectrum.  Also,  the  crosscorrelation, 
quadrature  and  co-spectra,  phase  difference,  coherence- 
squared  and  linear  regression  analysis  were  calculated 
between  specific  pairs  of  channels  (see  Table  7).  Examples 
of  cross-spectral  analysis  results  for  runs  two  and  12  are 
shown  in  Pigs.  12-15.  Additional  results  are  given  in 
Appendix  A.  Thermistor  2  on  the  horizontal  array,  TMH2 , 
was  selected  for  temperature  comparisons  with  the  sound 
records  because  this  thermistor  was  most  central  to  the 
sound  propagation  path  during  the  runs. 

1.  Basic  Statistics 

Representative  basic  statistics  of  temperature, 
waves,  sound  amplitude  and  sound  phase  are  shown  in  Table 
8.  By  examining  the  order  of  magnitude  of  these  statistics 
in  relation  to  the  depth  of  the  run  and  the  proximity  to 
the  thermocline,  certain  qualitative  relationships  between 
the  parameters  can  be  seen: 

1)  The  magnitude  of  the  temperature  variance  was 
determined  primarily  by  proximity  to  the  layer  and  its  mean 


TABLE  7 

CHANNELS  COMPARED  IN  TEMPORAL  STOCHASTIC  ANALYSES 


Abscissas  For  Regression  Analyses 
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FIGURE  14.  Run  2  SAMI  vs  TMH2  Spectra 


58 


Time  Series  Statistical  Parameters 
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period  was  controlled  by  surface  wave  swell  until  the  layer 
was  very  close  to  the  thermocline. 

2)  Sound  amplitude  variance  was  controlled  by  proximity 
to  the  thermocline  (temperature  gradient)  when  the  sonar 
array  was  horizontal.  The  movement  of  the  thermocline 

with  respect  to  the  sound  path  must  therefore  be  the  primary 
influence  on  sound  amplitudes  in  this  configuration. 

3)  The  amplitude  variance  when  the  array  was  nearly 
vertical  was  much  smaller  than  when  it  was  horizontal. 

This  indicated  the  effect  of  turbulence  proximity  to  the 
hydrophone  was  smaller  than  the  effect  of  temperature 
gradient  motion  over  the  sound  path. 

4)  Since  the  mean  period  of  the  sound  fluctuations  was 
much  smaller  than  the  mean  period  of  the  swell  at  depth, 
turbulent  interactions  related  to  the  higher  frequency 
internal  wave  motion  of  the  thermocline  were  emphasized 

in  the  sound  measurements. 

5)  When  the  sound  array  was  vertically  oriented  vice 
horizontal,  much  larger  phase  variations  occurred  even 
though  amplitude  variations  were  smaller.  This  resulted 
from  the  large  temperature  variations  over  thepath  length 
caused  by  waves  pumping  the  thermal  gradient  up  and  down. 

This  also  showed  that  the  temperature  variation  effects 
on  sound  amplitude  and  phase  were  in  quadrature. 

6)  The  order  of  the  mean  sound  phase  variances  indicates 
that  phase  was  primarily  influenced  by  surface  waves  when  the 
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array  was  horizontal.  Since  spectral  analysis  showed  that 
the  wave  coherence  with  temperature  was  very  low  at  the 
surface  wave  frequencies,  this  effect  may  be  bubble  related. 

2.  Spectral  Analysis 

Cross-spectral  analysis  of  the  temperature,  waves 
and  sound  records  for  the  horizontal  sonar  runs  one  and 
two  showed  that  the  sound  amplitude  changes  (Figs.  13  and 
It)  were  in  phase  with  the  waves  and  180°  out  of  phase  with 
temperature  in  the  regions  of  highest  coherence.  This  was 
expected  In  the  stable  density  structure  encountered. 
Additionally,  the  coherence  squared  function  (simply  re¬ 
ferred  to  as  coherence)  showed  that  only  the  long-period, 
prominent  swell  was  important  to  the  sound  amplitude  modu¬ 
lation  because  the  higher  frequency  wave  particle  velocities 
were  much  more  rapidly  attenuated  with  depth. 

Sound  phase  (Fig.  15)  was  found  to  vm\y  directly 
with  temperature  at  the  internal  wave  and  surface  wave 
frequencies  although  coherence  at  the  surface  wave  frequen¬ 
cies  was  much  smaller  than  at  the  internal  wave  frequencies. 
This  indicated  that  the  turbulence  associated  with  th-  layer 
was  more  important  than  bubbles  or  surface  wave  perturbation 
of  the  thermocline  in  influencing  sound  speed  over  the 
horizontal  sound  path.  A  comparison  of  the  coherences  of 
Figs.  12  and  13  shows  that  the  temperature  effect  is  at 
least  as  important  as  bubbles  even  at  the  surface  wave 
frequencies. 
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Similar  cross-spectral  analyses  with  the  sonar  path 
angles  25°  from  the  vertical  (run  12)  yielded  markedly 
different  results.  For  the  sound  path  entirely  above  the 
layer  (SAMI),  coherence  with  temperature  was  small  over  all 
frequencies.  Small  coherence  peaks  occurred  at  0.15Hz  and 
0.7Hz  which  indicated  that  the  vertical  motion  of  the  turbu¬ 
lent  blobs  forced  by  surface  swell  had  some  effect  In  the 
angled  sound  path. 

Sound  phase  coherence  with  temperature  was  largest 
at  the  Internal  wave  frequencies  but  smaller,  broad-banded 
coherence  at  the  surface  wave  frequencies  was  also  observed. 
This  confirmed  that  the  phase  was  being  controlled  by  the 
vertical  motion  of  the  thermocline. 

Most  of  the  sound  path  associated  with  SAM2  (2.3m) 
was  not  in  the  thermocline  during  the  majority  of  run  12. 
However  the  hydrophone  was  near  the  top  of  the  layer.  As  a 
result,  its  amplitude  was  in  phase  with  temperature  at  the 
internal  wave  frequencies  and  coherence  at  other  frequencies 
was  almost  zero.  This  verified  that  the  interference  effect 
of  the  proximity  of  the  thermocline-associated  turbulence 
to  the  receiver  was  theprimary  effect  on  SAM2  variability. 

Very  high  coherence  and  180°  phase  difference  between 
TMH2  and  SAM3  at  the  internal  wave  and  the  surface  wave 
frequencies  was  observed.  Since  the  acoustic  path  associated 
with  SAM3  generally  spanned  the  thermocline,  this  indicated 
that  the  variability  of  amplitude  across  the  entire 


63 


thermocline  was  controlled  by  refraction  and  diffraction 
similar  to  the  horizontal  case.  The  effect  of  diffraction 
was  smaller  than  for  the  horizontal  case,  because  the  len¬ 
ticular  shape  of  the  turbulent  blobs  present  broader  profiles 
to  the  vertical  sound  path.  The  reduced  diffraction  effects 
caused  the  variance  in  sound  amplitude  for  run  12  to  be  an 
order-of-magnitude  less  than  that  for  run  1. 

3.  Regression  Analysis 

p 

Sound  amplitude  and  phase  variances,  V  (SAM)  and 
V^SPH),  were  plotted  versus  TMH2  variances,  V2(T),  and 
least  squares  regression  analyses  were  applied  to  the  data. 
Linear  relationships  between  V2(SAM)  and  V^SPH)  versus 
V2 (T)  were  found  to  exist  for  runs  one  and  two  but  not  for 
run  12.  This  was  because  of  insufficient  variability  in 
the  portion  of  run  12  analyzed.  The  empirical  relationship 
below  was  developed  from  the  mean  regressions  calculated 
for  runs  one  and  two: 

V2(SAM)  -  (0.31»yb2/oC2)V2(T)e“L/x  (37) 

where  L  a  acoustic  path  length  and  x  =  the  distance  required 
to  reduce  V2(SAM)  by  1/e  (determined  empirically),  x  is 
near  the  size  of  the  Taylor  microscale. 

The  derivation  of  Equation  (37),  plots  of  the 
variances  over  the  run  segments  and  the  regression  plots 
are  given  in  Appendix  B. 
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C.  SPATIAL  ANALYSIS 


Temperature  spatial  correlation  and  structure  functions 
were  calculated  directly  from  ensemble  averages  of  26.67  min 
records  from  all  possible  combinations  of  thermistor  pairs 
along  the  vertical  thermistor  array.  Twenty-two  spatial 
lags  with  values  between  0  and  260  cm  were  thus  determined 
from  seven  thermistors.  The  value  of  C T  (see  Equation  33a) 
and  the  Taylor  microscale  as  well  as  the  Eulerlan  Integral 
scale  were  derived  from  the  above  functions  (Table  9). 

However,  since  large  variations  were  seen  in  these  functions 
the  validity  of  the  calculated  scales  was  questionable. 

In  an  attempt  to  remove  the  non-st ationary  effects 
resulting  from  internal  waves,  analyses  using  ensemble 
averages  of  short  segments  having  nearly  constant  turbulent 
intensity  were  conducted  for  runs  two  and  twelve.  The 
segments  were  chosen  so  that  examples  of  low,  medium  and  high 
turbulent  Intensity  could  be  examined  (Table  10) .  Because 
the  thermocline  steepness  was  highly  variable  along  the 
vertical  thermistor  array,  however,  only  the  medium  inten¬ 
sity  correlation  function  was  smooth  enough  to  provide  satis¬ 
factory  Eulerlan  integral  scales  and  the  structure  function 

of  all  samples  varied  widely.  Values  of  the  Taylor  microscale, 

2 

Eulerlan  integral  scale  and  CT  generated  from  hand-smoothed 
results  are  given  in  Table  11.  The  correlation  function  and 
structure  function  values  as  well  as  the  mean  variance  for 
the  vertical  array  thermistors  are  given  in  Appendix  C. 


TABLE  9 

Temperature  Spatial  Scales  (cm)  for  26.67  min  Ensembles 


R tn  * 

Taylor1 

Microscale 

Taylor^ 

Microscale 

Eulerlan 

Integral  Scale 

1 

15 

40 

36 

1.45x10"^ 

2 

178 

158 

175 

3.6xl0“3 

12 

60 

100 

36 

8.3xl0”3 

TABLE  10 

Relative  Turbulent  Intensity  for  Special  Spatial  Analyses 


Qualitative 

Intensity 

Enseirble 
Size  (min) 

TME  Mean 
Variance (°C)2 

Range 

(°c? 

Run 

Low 

4.67 

0.27X10-1 

±0.77xl0”2 

12 

Medium 

4.67 

4.24X10"1 

±4.95xl0”1 

2 

High 

3.73 

7. 79xlO~1 

±2.23xl0"1 

2 

TABLE  11 

Temperature 

Spatial  Scales  (cm)  for 

Special  Ensembles 

Qualitative 

Intensity 

Taylor1  Taylor2 

Microscale  Microscale 

Eulerian  «2 

Integral  Scale 

Low 

23 

6.8 

24 

3*3x10" 

Medium 

71 

40 

35 

4.4x10" 

High 

15 

20 

63 

5.4x10" 

Taylor1  ■  Taylor  microscale  from  spatial  correlation 
Taylor2  ■  Taylor  microscale  from  structure  function 


D.  APPLICATION  OP  MODELS 

While  it  Is  apparent  that  some  of  the  scale  estimates 
calculated  from  the  vertical  thermistor  array  may  be  erron¬ 
eous,  It  is  instructive  to  apply  the  theoretical  models 
for  the  sound  amplitude  variance  to  see  which  model  is 
most  applicable  in  the  non-stationary  case.  Table  12  lists 
the  equations  necessary  to  calculate  the  basic  parameters 
for  application  of  the  models.  For  the  purposes  of  this 
examination,  the  values  of  the  Richardson  number  (R^),  the 
Prandtl  number  (o^)  and  u/cQ  were  assumed  to  be  the  same 
as  reported  by  Neubert  [1970]  (see  Table  13  for  values). 

The  speed  of  sound  is  used  to  calculated  Wilson’s  Equation 
[Urick  1967]. 

Tables  13  and  14  list  the  resultant  parameters  of  each 
sample  for  which  scales  were  determined.  It  is  noteworthy 
that  the  values  of  the  refractive  index  (a)  calculated  above 
are  an  order  of  magnitude  above  the  near  surface  results 
reported  by  other  investigators  [Haley  1972,  Medwin  1973] 

(see  Table  15)  except  in  che  low  turbulent  intensity  cases. 
This  is  further  evidence  that  the  relative  influence  of  the 
turbulent  thermocline  in  sound  scattering  is  much  greater 
than  bubbles. 

If  the  order  of  magnitude  of  the  Taylor  microscale  is 

correct  (see  Table  13  for  comparative  values)  the  values  of 
2 

a  a  are  one  to  four  orders  of  magnitude  greater  than  previously 
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TABLE  12 

EQUATIONS  USED  FOR  MODEL  PARAMETER  EVALUATION 


Parameter 

Equation 

co 

c  “  1449  +  4.6T  -  0.55T2  +  0.17d 

0 

a 

a2  =  [4. 62V2(T)  +  0.055V(T)]/co2 

V 

Xg  ■  ai^  B  2.65a 

ng  -  Ag4(Rt/15)  -  2.07  x  loV 

c„ 

C  2  =  C  2/4T  2 
n  T  '  xo 

» 

lo 

l  *  =  0(2. 2/2. 4a)  »  0(0. 9a) 

0 

T  =  mean  temperature  (°C) 

d  B  depth  (m) 

V2(T)  =  variance  in  temperature  (°C2) 
a  =  Taylor  microscale  (cm) 

n  =  Kolmogorov  microscale  (cm) 

s 

*  equation  for  lQ  assumes  u/c0  Is  a  constant 
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TABLE  13 

Experiment  Parameter  Values  for  Model  Tests 


Parameter 

c0(m/sec) 


L  (cm) 
a  (cm) 


X  (cm) 
kQ  (cm"1) 


RUN  1 
1516 

4.5x10"^ 


13.5 


0.84 

1.45x10' 


RUN  2 

1498 

2.1x10' 


0.84 

1.05x10' 


RUN  12 

1515 

7.9x10' 


0.84 

2.4x10' 


Le  (cm) 


Xg  (cm) 

ne  (m) 


10.5 

6300 

3.0x10' 


21.4 

6300 

7.9x10' 


19.2 

6300 

3.7xlO"5 


—  m'i  jn'i  1 


TABLE  14 


Short  Sample  Parameter  Values  for 

T.nw  Medium 

Model  Tests 

High 

Parameter 

MW  ** 

Intensity 

Intensity 

Intensity 

c_  (m/sec) 

1512 

1508 

1500 

0 

a 

5.0xl0~4 

2.0xl0”3 

2.7xl0“3 

L  (cm) 

.  200 

200 

200 

a  (cm) 

23 

71 

15 

o 

* 

21 

64 

13 

X  (cm) 

7.5 

7.5 

7.5 

k  (cm”1) 

0.84 

0.84 

0.84 

0 

c 

-4 

1.05x10 

I.l6xl0”4 

1.29x10” 

un 

Le  (cm) 

24 

35 

63 

V 

°k 

7 

7 

7 

X„  (cm) 

61 

188 

40 

E 

n  (cm) 

2.8 

8.5 

1.8 

6 

R 

6300 

6300 

6300 

"t 

a2  a  (cm) 

5. 8xl0~6 

2. 8xl0”4 

l.lxlO-4 
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reported.  This  result  has  a  direct  relationship  (in  the 
Born  model)  to  the  variance  of  sound  amplitude. 

The  variances  of  sound  amplitude  calculated  by  each 
model  previously  discussed  (Geometric,  Born,  Rytov,  Debye) 
are  presented  in  Tables  16  -  19.  In  general,  the  best  fit 
to  the  measured  sound  variances  were  calculated  from  the 
Born  model  using  the  Krasil’nikov-Obhukov  formulation.  A 
surprisingly  good  fit  to  the  true  variance  was  also  obtained 
from  the  geometric  model  in  the  mid-intensity  range.  Because 
of  the  above,  the  stimates  of  the  Taylor  microscale  seem 
reasonable. 

2 

Since  the  accuracy  of  the  scales,  and  L&,  the  Eulerian 
integral  scale,  is  in  question,  no  firm  statement  can  be 
made  about  the  adequacy  of  the  Rytov  and  Debye  models. 

However,  because  L  is  thought  to  be  correct  for  the  mid- 

C 

intensity  ensemble,  the  Debye  model  appears  incorrect. 
Presumably,  calculation  of  the  structure  function  and  corre¬ 
lation  function  from  the  thermistors  in  the  horizontal  array 
over  stationary  ensembles  would  resolve  this  question. 
Unfortunately,  this  was  not  accomplished  in  time  for 
inclusion  in  this  thesis. 


TABLE  16 

Geometric  Model  Experimental  Results 


V2  ■  4/15  /if  O2 

(L/a)3  .  i 

7.5  «  l02fr 

(cm) 

Ensemble 

Model  V2 

Observed  V2(yb) 

2  lQ2/l  (cm) 

Run  1 

4.8xl0"4 

1.12xl0“4 

0.91 

Run  2 

'  3xl0”6 

1.45xl0”3 

128 

Run  12 

l.lxlCf5 

4.05xl0~5 

146 

Low  Intensity 

7.8x10”5 

7.1xl0”6 

2.2 

Mid  Intensity 

4.2xlO"5 

9.5xlO“5 

21 

High  Intensity 

8.2xl0“3 

3. OxlO”4 

0.85 

TABLE  17 

Born  Model  Experimental  Results 
[Krasil'nikov  and  Obhukov] 


V2  -  k02a2aL(l  - 

1/Dtan"1  D) 

j  D  *  irL/kQa 

Ensemble 

Model  V2 

2 

Observed  V 

_ D 

Run  1 

4.3xl0~4 

1. 12xl0”4 

50 

Run  2 

7.9xlO"4 

1.45xl0”3 

4.2 

Run  12 

3.74x10"5 

4.05xl0-5 

12.5 

Low  Intensity 

5.8x10‘6 

7.1xl0-6 

32.5 

Mid  Intensity 

2 . 8xlO”4 

9.5xl0"5 

10.5 

High  Intensity 

_li 

1.1x10 

3.0xl0"4 

50 

D  >>  1 

Wave  Region 

TABLE  18 


Rytov  Model  Experimental  Results 
V2  -  161  Cn2  L  11/6  -7/6  ;  7.5  « 


Ensemble 

Model  V2 

2 

Observed  V 

sA^r 

Run  1 

5.2xl0“3 

1.12xl0“4 

5.5 

Run  2 

2. 8x10" 3 

1.45xlO“3 

150 

Run  12 

1.5xl0"2 

4.05xl0”5 

35 

Low  Intensity 

2. 8xl0“3 

7.1xl0“6 

9.9 

Mid  Intensity 

3.4x10“ 3 

9.5xlO“5 

44 

High  Intensity 

4. 2xl0“3 

3. 0x10“^ 

5.2 

TABLE  19 

Debye  Model  Experimental  Results 

V2  -  2c2k02LeL 

Ensemble 

Model  V2 

? 

Observed  V 

Run  1 

2.1xl0“3 

1.12X10"4 

Run  2 

0.22 

1.45xlO“3 

Run  12 

6.3xlO”3 

4.05xl0‘5 

Low  Intensity 

1.7xl0"3 

7.1xl0"6 

Mid  Intensity 

4xl0“2 

9.5xlO“5 

High  Intensity 

0.13 

3.0x10”^ 
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V.  CONCLUSIONS 


The  theories  of  turbulence  which  have  been  developed 
for  the  ocean  require  some  degree  of  stationarity  in  the 
medium  before  they  may  be  applied.  This  paper  has  examined 
four  basic  stochastic  models  of  acoustic  wave  propagation 
under  the  most  severe  conditions  of  non-stationarity,  e.g. 
near  a  sharp  thermocline  which  is  being  oscillated  by  both 
large  surface  swell  and  by  Internal  waves. 

The  experimental  conditions  specified  that  the  solution 
form  would  be  in  the  wave  region  for  which  the  Born  model 
would  be  most  accurate.  The  solution  of  the  Krasil’nikov- 
Obhukov  form  of  the  Born  model  from  directly  measured, 
ocean  parameters  gave  reasonable  results. 

The  solution  of  the  Rytov  and  Debye  models  for  the  same 

conditions  yielded  results  two  orders  of  magnitude  higher 

than  was  actually  measured.  Because  evaluation  of  scales 

was  difficult,  this  was  expected  even  though  both  models 

are  more  general  than  the  Born  model  and,  in  fact,  contain 

its  solution  in  a  subset.  This  also  showed  the  determination 

of  the  Taylor  microscale  was  more  reliable  than  the  deter- 

2 

mination  of  the  structure  constant,  C,p  ,  or  the  Eulerian 
integral  scale  for  measurements  across  the  thermocline. 

An  attempt  was  made  to  apply  the  correlation  function 
and  structure  function  piecewise  to  small  sections  of  data 
which  were  thought  to  be  relatively  stationary.  For 
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application  across  the  thermocline,  this  method  did  not  work 
well  because  selection  of  such  sections  Is  very  difficult. 

The  results  of  the  time  series  analyses  showed  that  the 
losses  over  the  acoustic  path  were  greatest  when  the  array 
was  horizontal  near  the  thermocline  because  the  diffractive 
losses  were  an  order  of  magnitude  greater  than  those  for 
the  nearly  vertical  case.  Interference  effects  due  to  the 
proximity  of  turbulence  to  the  receiving  hydrophone  accounted 
for  a  0.14?  signal  loss  and  scintillations  due  to  surface 
swell  accounted  for  0.07?  signal  modulation.  The  effects 
of  bubbles  were  small  because  the  measurements  were  rela¬ 
tively  deep  and  the  acoustic  frequency  was  much  lower  than 
the  peak  bubble  resonant  frequency.  Measured  loss  over  the 
2m  path  was  on  the  order  of  5?  and  phase  shifts  of  up  to 
30°  occurred  for  these  worst-case  conditions. 

Current  sound  propagation  models  based  upon  empirical 
data  have  undoubtedly  included  these  effects  in  surface 
scattering  or  reverberation.  But  the  roughness  description 
for  the  sea  surface  used  by  sonar  models  is  based  upon  sea 
state.  The  losses  here  described  are,  instead,  related  to 
the  swell  wave-induced  water  particle  velocities  and  internal 


waves. 


wmmmmm 


APPENDIX  A 

Spectral  Analysis  Plots 

This  appendix  includes  examples  of  cross-spectral 
analysis  plots  so  that  qualitative  comparisons  between 
sound,  temperature  and  waves  during  the  experiment  can  be 
made.  Included  are  spectra  for  1)  sound  amplitude  versus 
temperature,  2)  sound  phase  versus  temperature,  3)  wave- 
related  effects  and  *0  comparison  of  one  sound  channel 
with  another  for  the  same  run.  The  plots  included  are 
listed  below. 

Sound  Amplitude  Versus  Temperature 


Run  # 

1 

1 

1 

2 

2 

2 

12 

12 

12 


Channels  Compared 

SAMI  -  TMH2 
SAM2  -  TMH2 
SAM3  -  TMH2 

SAMI  -  TMH2  (see  text) 
SAM2  -  TMH2 
SAM3  -  TMH2 

SAMI  -  TMH2 
SAM2  -  TMH2 
SAM3  -  TMH2 


Sound  Phase  Versus  Temperature 


Run  # 


2 

2 


12 


Channels  Compared 
SPH1  -  TMH2 

SPH1  -  TMH2  (see  text) 
SPH2  -  TMH2 

SPH1  -  TMH2 
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APPENDIX  B 

Analyses  of  Variance  Regressions 

Variances  of  56  sec  samples  of  data  which  included  1400 
data  points  each  were  calculated  along  26.67  min  segments 
of  certain  channels  of  runs  one,  two  and  12.  Regression 
analyses  for  sound  amplitude  and  phase  versus  temperature 
were  computed  using  these  variances.  The  resultant  plots 
of  the  variances  with  the  regression  curves  are  shown  after 
Table  B5.  A  summary  of  the  regression  equations  is  given 
in  Table  Bl.  The  following  mean  variance  equations  were 
computed  from  the  average  of  the  regression  equations  listed 
for  runs  one  and  two: 


V2(SAM1)  =3.2 

x  10"3  V2(T) 

(Bl) 

V2(SAM2)  =  1.6 

x  10‘3  V2(T) 

(B2) 

V2(SAM3)  =  4.1 

_1i  p 

x  10  v  (T) 

(B3) 

If  V2(T)  is  eliminated  from  the  above  equations  (stationary 
assumption),  then  because  SAMI  and  SAM2  are  30  cm  apart  and 
60  cm  separates  SAM2  and  SAM3: 

V2(SAM)  cC  (i)"L/3°  cm  (B4) 

where  L  =  sound  path  length. 
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This  equation  can  be  checked  by  comparing  the  sound 
amplitude  variances  for  specific  56  sec  samples  (where  the 
sample  Is  stationary)  and  similarly  deriving  the  relation¬ 
ship 

V2(SAM)  oc  (|)"L/x  (B5) 

where  x  is  the  distance  required  to  reduce  V2 (SAM)  by  0.5. 
This  is  done  in  Tables  B3  -  B5.  The  following  equations 
result : 

9  v-L/34.5 

Run  1:  V*(SAM)<*  ±  cm  (B6) 

i 

Run  2:  V2(SAM)<x  cm  (B7) 

Conversion  from  base  j  to  base  e  and  placing  the  path 
length  reference  at  lm  vice  2m  yields  (from  equations  Bl, 
B4,  and  B5) : 

V2(SAM)  =  (0.3^ub2/°c2)V2(T)e”L/x  (B8) 

where : 

L  =  acoustic  path  length  in  meters, 
x  ®  temperature  spatial  scale  for  which 


V2(SAM)  -  V2  (SAM)x>10/o 
o 

V  (T)  ■  temperature  variance 

2 

V  (SAM)  *  variance  of  sound  amplitude. 

And  the  spatial  scales  in  (B4),  (B6)  and  (B7)  become  ^cm, 

50cm  and  6lcm,  respectively. 

2  2 
Since  V  (SAM)  is  known  to  be  a  function  of  0  »  the 

variance  of  the  index  of  refraction,  and  a  “  2.1  x  10 

for  run  two  (which  had  the  best  regression  fits): 

V2(SAM)  =  7.7  x  10^  a 2  V2(T)  e“L/x  (B9) 

Assuming  that  x  <<  L,  the  following  approximation  can  be 
made: 


x  „  x 
S+L  "  L 


(BIO) 


Substituting  Equation  B ( 10 )  into  Equation  B(9)  and  multi- 

p 

plying  by  L  to  compensate  for  the  spherical  spreading  loss 
yields: 


V2(SAM)  =  7.7  x  101*  a2  V2(T)  xL 


(Bll) 


which  is  similar  to  the  exponential  form  of  the  Bom 
approximation: 

V2  -  k  2  a2  aL  (B12) 

o 

if  kQ  is  similarly  included  in  Equation  (Bll): 

V2(SAM)  -  1.1  x  10“5  kQ2  a2  V2(T)  xL.  (B13) 

The  plots  following  Table  B5  are  included  to  allow 
qualitative  comparisons  of  conditions  between  runs  and  of 
the  variances  between  channels  during  the  same  run.  For 
each  run  analyzed  the  variance  versus  time  plots  are  given 
followed  by  regression  plots  between  various  channels. 
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TABLE  Bl.  Regression  Equations  for  Sound  Versus  Temperature 
Variances 

AMPLITUDE 


Channel 


SAMI 


Run# 


SAM2 


SAM3 


Equation 

V2(SAM1)  -  [B.W^T)]  x  10"3 

V2(SAM1)  “  C3.0V2(T)]  x  10"3 

V2 (SAMI)  '  [4  -  2V2(T) ]  x  10“4 

V2 (SAM2)  =  [2.0V2(T)3  x  10”3 

V2(SAM2)  ■  [1.2V2(T)]  x  10"3 

V2(SAM2)  *  [0.7  +  1.5V2(T)]  x  lO-5 

V2 (SAM3)  *  [4. 8V2(T) ]  x  10’4 
V2(SAM3)  -  C3.3V2 (T) D  x  10-4 
V2(SAM3)  3  [0.2  +  4V2(T)]  x  lO"4 
PHASE 


Channel 

SPH1 


SPH2 


Run#  Equatlon 

1  V2(SPH1)  =  1  +  *UV2  (T) 

2  V2 (SPHl)  =  0.2  +  0.7V2 (T) 

12  V2(SPH1)  -  3.7  +  77V2 (T) 

1  V2(SPH2)  =  2  +  8V2(T)» 

2  V2(SPH2)  -  2  +  22V2 (T) 

12  V2(SPH2)  =  2  -  25V2 (T) 


•From  20  minute  record. 


TABLE  B2.  Run  1  Variances 


variance 

26  2/3  min 

56  sec 

56  sec 

56  sec 

v/^C2) 

2.2  x  10"2 

16.1  x  10"2 

3.2  x  10"2 

0.43  x  10* 

VSAMl(yb2) 

11.0  x  10’5 

52.0  x  10-5 

12.0  x  lO"5 

18  x  lO”5 

vsL>fc2) 

7.6  x  10"5 

39.0  x  10“5 

7.2  x  10“5 

13  x  10“5 

VSAM3(yb2) 

2.5  x  10-5 

9  x  IQ"5 

3.3  x  lO"5 

3.2  x  10“5 

TABLE  B3.  Sound  Amplitude  Distance  Loss  Factors  (Run  1) 


sanple 

distance 

30  cm 
60  cm 
90  cm 

quantity 


26  2/3  min 

1. 45 
3.04 
4.40 


Rm  1-1 
56  sec 

1.33 

4.33 
5.78 


standard 

deviation 


Run  1-2  Run  1-3 

56  sec  56  sec 

1.67  1.38 

2.18  4.06 

3.64  5.63 

lormalized  accepted 


distance 

factors 

in  loss 
factors 

factors 

factors 

30  cm 

1.46 

0.18 

1.00 

1.15 

60  .cm 

3.52 

1.17 

2. 4l 

2.30 

90  cm 

5.02 

1.19 

3.44 

3.45 

•not  included  in  mean  loss  factors 


TABLE  B4.  Run  2  Variances 


Sanple 

An  2 

Rtn  2-9 

An  2-26 

An  2-20 

variance 

26  2/3  min 

56  sec 

56  sec 

56  sec 

VjW) 

46.2  x  10"2 

98.1  x  10“2 

143.4  x  10“2 

121.4  x  10-2 

VSAML(yl>2) 

140  x  10“5 

311  x  lO"5 

422  x  lO"5 

420  x  lO”5 

VsLb<*2> 

56.0  x  10-5 

122  x  10“5 

146  x  10"5 

178  x  10"5 

CM 

35 

s> 

17.0  x  10*5 

30.5  r  10-5 

40.3  x  lO"5 

460  x  lO"5 

TABLE  B5. 

Sound  Amplitude 

Loss  Factors  with  Distance  (Run  2) 

sample 

Sun  2* 

Rum  2-9 

An  2-26 

is* 

distance 

26  2/3  min 

56  sec 

56  sec 

56  sec 

30  cm 

2.5 

2.5 

2.9 

2.4 

60  cm 

3.29 

4.0 

3.62 

3.9 

90  cm 

8.24 

10.2 

10.47 

9.1 

quantity 

distance 

mean 

loss 

factors 

standard 

deviation 

In  loss 
factors 

normalized 

loss 

factors 

accepted 

loss 

factors 

30  cm 

2.6 

0.3 

1.0 

1.4 

60  cm 

3.8 

0.20 

1.5 

2.8 

90  cm 

9.9 

0.73 

3.8 

4.2 

•not  Included  In  mean  loss  factors. 


ft  UKj  l 


x-scmr.2. 001  *02  UNITS  INCH,  --  Timb  u«o 
y-5>CfiLF--2.00f.-02  UNITS  INCH.  * VW)  C*c») 

URRIFINCL  TIT 


X-St;nu  •••J.OOr+02  UNITS  inch. 
y-SCfilf  •i'OOE -0!)  UNITS  INCH. 

UflRIRHCt  FIT 


{.sbc) 


0J2 


CIO 


o:i  005  qc‘j 

RunL  I 

K-SCnL  E -2- 00f >02  UNITS  INCH.  =  Time  (sec) 
v'-Sffllf  ■  i  •  001. -OH  UNITS  INCH.  =  Vl(srt«a) 

URRIRNCE  FIT 


109 


Run!  I 


X-SCRLF.---2.00F  -02  UNITS  INCH. 
^'-Sr.RL  F'-l«  OOF -OS  UNITS  INCH. 

REGRESSION  PLOT 


111 


(6CX>) 

V1  (  SAMI  )  (VOLTS’1  ) 


"•^•Mssseiai 


I 

k 


Run)  I 


K- Sf.fiL f  2 •  OOF  -02  UNITS  INCH. 
V-SCRLF  ■-!- OOF  -*00  UNITS  INCH. 

RKCRKSSION  PLQT 


r  V1  (TMHX)  C  °Cl) 

-  V1  C  SPH  1 )  ^  decrees' ) 


11* 


ummmiUMiM 


mum**** 


nn'-i  i-iii^JiMliMiiiiiailiMfe^aa 


K-5CPHF.  -2.00F-02  UNITS  INCH. 

y-scmr  -5.  oor-02  units  inch. 


REGRESSION  PLOT 


-  V2(  TAAH  1  )  (°Cl) 

-  vxC  ^ 


115 


c  -- 


K -r/ni  f  ■ •  OOF <02  UNITS  INCH 
v-SffiiF'-U.OUf-  02  UNITS  MM- 

URRIRNCl  FIT 


r  TIME  C  ^Gc  ) 

=  vz(wrtve‘)  t»'w 


116 


fCC  see 


*■»»> 


RUNiX 


K-ST.RLF--2.00F-.02  UNITS  INCH  -  Tim .v 
v-3C.nU-2.lWf  01  UNITS  INCH  -M  (rMMa>  (  c  ' 

UfiR I RhCt  FIT  117 


Run/  X 


K-SCAIF  -?.00f->02  UNITS  INCH.  •  r.*«ecs«c'> 

Y-Sr.i-1L c  1  0Cf-O<J  UNITS  INCH  W«"-Ti  ) 

URR I fiNCL  FIT 


Ruw  X 


X-SCP.IF-  2. 00!' +02  UNITS  INCH.  -  time  Uec) 
y-SCPLF'--i. OOF-04  UNITS  INCH.  *  l 

URR1RNCE  FIT 


119 


R«k)  X 


X-SC.RU  " 2- POE  +02  UNITS  INCH. 
y-SCRl !'  :S. 02F-06  UNITS  inch. 

URRIRNCF  FIT 


120 


s-srmr- 

y- SCALE 


'  i.lilf+02  UNITS  INCH-*  TiMe(»ec)  . 
S.OOt-O.l  UNII3  INC.H.  =  Vx(sphi)  U^re.i  ) 


URR1RHCE  FIT 


121 


<u»y,  t; 


X-SCfiLfNf-POr+02  UNITS  INCH.  -  TIMe  l iecT 
Y-SCP>U>5.00t>00  UNITS  INCH.*  VHspiu) 

URRIRNCF  FIT  122 


1211 


R  U  Nl  X 


3CRL  T---2- OOt'-Ol  UNITS  INCH* 

mr-^.oor-os  units  inch. 

■GRES5I0N  PLOT 


*  va(  r^Hi)  (°cz) 

-  v2(  SA/v\3)  (vOUTSj 


ri-Sr.Riri’.Ot1?  0.!  UNITS  INCH*  Va(rMH2)  t°c') 

Y-3r.ni  r  ■R.GOP  U!  UNITS  INCH- *  VMSPHi)  (dejfees  ) 

REGRESSION  PLOT 

126 


STS  OTO 


§ 


K-mr--2.bOF.-oi  units  inch.  =  (*0 

y-Sr.Rl r--5.00F+00  UNITS  INCH.  =  (deyees) 

REGRESSION  PLOT 


127 


RUtJ  12- 


K-SCRLJ  ■?. OOf -*02  UNITS  INCH.  =  PmeUcc) 
Y-SCAl.r-.-5.0OE -02  UNITS  INCH.  =  V’-(wame')(, 

URRIRNCt'  FIT 


K-Sr.HL F  -2- OOf  +02  UNITS  INCH.  »  T./we  Csec) 

V-5C-RL F -•-&* OOf  -06  UN3T3  INCH.  *  Vllsrt«3)  (volts1) 

URRIflNCF  FIT 


X-SCFIl  E  ■-2.00E+02  UNITS  INCH.  -  riwe  (ssc  )  t 

w_S[.fliE  1.  OUf.+Ol  UNITS  INCH  -  -  vMsphO  (‘JeTr**s  ) 

UflR I RHCE  FIT 

13*1 


II,  -u > jau^ijpa^g^iiWBip^BPIPPgllBBgiWBgBBgpSBIIPWWPIIB^P^ 


tfMMM*/'  »<**> 


Run)  IX 

X-SRRl f. - 1  •  OOF  -0.1  UNITS  INCH.  =  n^ItmHz)  (°cl) 
V-Sr,Rlf  -5.001-06  UNITS  INCH-  *  -vM'-rtMi  )  <  v*utC,  ) 

REGRESSION  PLOT  ■ 


RukJ  IX 


K-SC.PiL E  -J.-  OOf.-Ol  UNITS  INCH.  *  V*(  mtu')  C°cx) 
y-SC.nLT-5.0D'  -0?  UNITS  INCH.  -  (volts*) 

REGRESSION  PLOT 


Rum  ix 

K-5C.ni  f  ■  1  DOF  -0.1  UNITS  INCH.  *  V2  C  C*c*) 
V-SCRLF'  S-OOF -OR  UNITS  INCH.  =  v2  (so*.  3)  (vouts*) 

RtGRKSS I  ON  PLOT 


rum  ix 


K-SC.RLF  -  l.  OOr-OJ  UNITS  INCH.  =  v*  Ctmhi')  (°c‘) 
V-SCFiLf  7-1 . 00!  +01  UNITS  INCH.*  V'CSPHl)  (press') 

REGRESSION  PLOT 


139 


Run)  IX 


X-SCALEV  J •  OOF-Ol  UNITS  INCH.  *  vMtmhi)  (  "c1) 
-'-SHALE' -  Ij. 00f-02  UN115  INCH.  -  v1  (  sPHa)  (decree 

REGRESSION  PLOT  lhn 


am* 


Run) 


K-Smif  1- OOf  -  G1  UNITS  INCH.  =  ^  \ 

^-Sr.FlL r -=s. OOf-02  UNITS  INCH.  =  V  (TMHP  (  «-  ) 

REGRESSION  PLOT  lU 


APPENDIX  C 


Spatial  Correlation  and  Structure  Functions 

\ 

The  value^  of  the  spatial  correlation  function  and 

i 

structure  function  calculated  during  spatial  stochastic 
analyses  are  tabulated  in  Tables  Cl  -  C6.  The  functions 
were  calculated 'along  the  vertical  thermistor  array. 

The  variance  of  Ihe  temperature  at  each  thermistor  over 
the  ensemble  is  a!(,so  given. 

The  Tables  include: 


Table  # 

Ensemble 

Cl 

.  Run  1,  26.67  min 

C2 

Run  2,  26.67  min 

C3 

* 

Run  3>  26.67  min 

Ck 

;  Low  turbulence  intensity 

C5 

Middle  turbulence  intensity 

C6 

High  turbulence  intensity 

TABLE  Cl.  Run  1,  Vertical  Spatial  Correlation  and 
Structure  Function. 


Thermistor  No. 

Variance 

(°C2) 

1 

0.5999 

2 

1*0970 

3 

O.0540 

4 

0.3090 

5 

0.2722 

6 

0.1359 

7 

O.0747 

Correlation 

Structure 

Lag 

Function 

Function 

— 

1.0000 

0.0 

0.0 

0.9635 

0.2081 

5.0000 

0.3745 

0.5841 

15.0000 

0.3654 

0.4673 

20.0000 

0.5507 

0.3978 

25.0000 

2.6243 

0.9860 

25.0000 

0.5208 

0.5926 

?o.oooo 

0.1736 

0.3503 

45.0000 

0.2716 

0.2302 

55.0000 

0.2668 

0.7407 

80.0000 

0.2430 

0.9232 

85.0000 

0.0721 

0.3287 

100.0000 

0.5209 

3.1869 

135.0000 

0.9274 

3.3103 

150.0000 

0.9724 

3.0852 

155.0000 

0.3622 

1.6778 

160.0000 

0.6725 

1.9936 

175.0000 

0.6874 

1.8289 

180.0000 

0.5094 

3.5094 

180.0000 

0.3640 

1.9750 

205.0000 

0.1944 

3.9140 

235.0000 

0.1514 

2.1753 

260.0000 
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TABLE  C2.  Run  2.  Vertical  Spatial  Correlation  and 

Structure 

Function 

Thermistor  No. 

Variance 

c\» 

0 

0 

1 

0.2154 

2 

0.1750 

3 

0.3544 

4 

0.5783 

5 

0.3910 

6 

0.4338 

7 

0.5170 

• 

Correlation 

Structure 

Lag 

Function 

Function 

1.0000 

0.0 

0.0000 

1.0592 

0.3607 

5.0000 

1.0090 

0.3676 

15.0000 

0.9452 

0.0160 

20.0000 

1.0354 

0.0388 

25.0000 

1.0891 

0.4103 

30.0000 

0.9706 

0.0749 

45.0000 

1.1503 

0.1313 

55.0000 

1.0532 

0.2159 

80.0000 

1.0844 

0.6326 

85.0000 

0.9813 

0.2662 

100.0000 

0.5463 

0.2614 

135.0000 

0.6164 

0.6939 

150.0000 

0.5559 

0.3357 

155.0000 

0.5010 

0.4553 

160.0000 

0.5551 

0.9197 

175.0000 

0.4985 

0.5539 

180.0000 

0.5580 

0.4413 

180.0000 

0.4878 

0.6848 

205.0000 

0.5458 

0.6783 

235.0000 

0.4557 

0.9619 

260.0000 

TABLE 


C3.  Run  12,  Vertical  Spatial  Correlation  and 
Structure  Function 

p 

Thermistor  No.  Variance  (°C  ) 

1  0.2792 

2  0.2433 

3  0.1690 

4  0.1444 

5  0.1275 

6  0.0743 

7  0.0083 

Correlation  Structure  Lag 

Function  Function 

1.0000 

0.0 

0.0 

1.0558 

0.0211 

5.0000 

1.1642 

0.1301 

15.0000 

1.0859 

0.1540 

20.0000 

0.7096 

0.1640 

25.0000 

1.5349 

1.0283 

25.0000 

0.7463 

0.2231 

30.0000 

0.7552 

0.3984 

45.0000 

0.1447 

0.3588 

55.0000 

0.1852 

0.6946 

80.0000 

0.1969 

0.8039 

85.0000 

0.1995 

0.9920 

100.0000 

0.8057 

1.6228 

135.0000 

0.6778 

1.6854 

150.0000 

0.6295 

1.6494 

155.0000 

0.6143 

2.2873 

160.0000 

0.5148 

2.2931 

175-0000 

0.4773 

2.2355 

180.0000 

0.4121 

I.6676 

180.0000 

0.2998 

2.1737 

205.0000 

0.1487 

1.6762 

235.0000 

0.1214 

2.0123 

260.0000 
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TABLE  C4. 


Low-Intensity,  Vertical  Spatial  Correlation 
and  Structure  Function 


Thermistor  No. 

Variance  (°C2) 

0.0666 

X 

p 

0.0774 

C 

0.0062 

5 

ll 

0.0038 

■t 

c 

0.0028 

0 

C 

0.0114 

u 

7 

0.0014 

Correlation 

Function 

Structure 

Function 

Lag 

1.0000 

0.0 

0.0 

1.0785 

0.2173 

5.0000 

1.3019 

0.9860 

15.0000 

1.0129 

1.1896 

20.0000 

0.2793 

4.5337 

25.0000 

9.4721 

32.6599 

25.0000 

0.2274 

5'.  0118 

30.0000 

0.0791 

6.1497 

45.0000 

0.1556 

4.2863 

55.0000 

0.4158 

0.6737 

80.0000 

0.4850 

0.9095 

85-0000 

0.6269 

1.4669 

100.0000 

3.1748 

23.5975 

135.0000 

1.6307 

25.8445 

150.0000 

1.5337 

25.6640 

155.0000 

2.9846 

20.1207 

160.0000 

1.5724 

22.1038 

175.0000 

1.1432 

22.5880 

180.0000 

0.7332 

30.3575 

180.0000 

-1.6798 

31.3264 

205.0000 

0.5584 

27.1199 

235.0000 

1.1856 

22.0083 

260.0000 
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TABLE  C5.  Mid-Intensity,  Vertical  Spatial  Correlation 
and  Structure  Function 


Thermistor  No. 

1 

2 

3 

4 

5 

6 
7 


Variance  (°C2! 

0.3603  x  10“5 
0.3726  x  10”? 
0.8536  x  10”1 
0.1329 

0.1011  . 
0.8996  x  10“ ; 
0.3400  x  10”1 


Correlation 

Function 

Structure 

Function 

Lag 

1.0000 

0.0 

0.0 

1.0345 

0.2462 

5.0000 

0.9265 

0.3068 

15.0000 

0.8702 

0.1042 

20.0000 

0.8345 

0.2212 

25.0000 

0.0002 

0.0034 

25.0000 

0.8742 

0.4570 

30.0000 

0.6962 

0.3425 

45.0000 

0.2549 

0.7168 

55.0000 

0.2131 

0.9103 

80.0000 

0.2374 

1.1769 

85.0000 

0.1709 

0.8393 

100.0000 

0.0130 

0.8224 

135.0000 

0.0205 

.1.2779 

150.0000 

0.0134 

0.9769 

155.0000 

-0.0008 

0.8464 

160.0000 

-0.0005 

1.3162 

175.0000 

-0.0008 

1.0017 

180.0000 

0.0149 

0.8641 

180.0000 

-0.0004 

0.8911 

205.0000 

0.0154 

0.3093 

235.0000 

0.0005 

0.3356 

260.0000 
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TABLE  C6.  High-Intensity,  Vertical  Spatial  Correlation 
and  Structure  Function 


Thermistor  No. 


Variance  (°C*) 

0.1800  x  10"5 
0.2462  x  10"| 
0.4578  x  10”? 
0.1231  x  10~i 
0.1226  x  10’; 
0.5298  x  10”f 
0.2655  x  10"1 


Correlation 

Function 

Structure 

Function 

Lag 

1.0000 

0.0 

0.0 

0.6971 

9.6513 

5.0000 

0.3813 

9.6565 

15.0000 

0.5125 

0.3485 

*20.0000 

1.7985 

1.7245 

25.0000 

0.0147 

0.1730 

25.0000 

1.1782 

12.0107 

30.0000 

0.8521 

2.9908 

45.0000 

8.3755 

9.2314 

55.0000 

3.7301 

15.2006 

80.0000 

2.9532 

25.7984 

85.0000 

1.9051 

18.2240 

100.0000 

0.2184 

0.1375 

135.0000 

0.2522 

9.7420 

150.0000 

0.2551 

0.6906 

155.0000 

0.0156 

0.3438 

160.0000 

0.0156 

10.0158 

175.0000 

0.0170 

0.9675 

180.0000 

0.3564 

3.8096 

180,0000 

0.0179 

4.2871 

205.0000 

1.0541 

19.7533 

235.0000 

0.0543 

21.5538 

260.0000 
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